Formulating genome-scale kinetic models in the post-genome era.

The biological community is now awash in high-throughput data sets and is grappling with the challenge of integrating disparate data sets. Such integration has taken the form of statistical analysis of large data sets, or through the bottom-up reconstruction of reaction networks. While progress has been made with statistical and structural methods, large-scale systems have remained refractory to dynamic model building by traditional approaches. The availability of annotated genomes enabled the reconstruction of genome-scale networks, and now the availability of high-throughput metabolomic and fluxomic data along with thermodynamic information opens the possibility to build genome-scale kinetic models. We describe here a framework for building and analyzing such models. The mathematical analysis challenges are reflected in four foundational properties, (i) the decomposition of the Jacobian matrix into chemical, kinetic and thermodynamic information, (ii) the structural similarity between the stoichiometric matrix and the transpose of the gradient matrix, (iii) the duality transformations enabling either fluxes or concentrations to serve as the independent variables and (iv) the timescale hierarchy in biological networks. Recognition and appreciation of these properties highlight notable and challenging new in silico analysis issues

Using in silico models to simulate dual perturbation experiments: procedure development and interpretation of outcomes.

BackgroundA growing number of realistic in silico models of metabolic functions are being formulated and can serve as 'dry lab' platforms to prototype and simulate experiments before they are performed. For example, dual perturbation experiments that vary both genetic and environmental parameters can readily be simulated in silico. Genetic and environmental perturbations were applied to a cell-scale model of the human erythrocyte and subsequently investigated.ResultsThe resulting steady state fluxes and concentrations, as well as dynamic responses to the perturbations were analyzed, yielding two important conclusions: 1) that transporters are informative about the internal states (fluxes and concentrations) of a cell and, 2) that genetic variations can disrupt the natural sequence of dynamic interactions between network components. The former arises from adjustments in energy and redox states, while the latter is a result of shifting time scales in aggregate pool formation of metabolites. These two concepts are illustrated for glucose-6 phosphate dehydrogenase (G6PD) and pyruvate kinase (PK) in the human red blood cell.ConclusionDual perturbation experiments in silico are much more informative for the characterization of functional states than single perturbations. Predictions from an experimentally validated cellular model of metabolism indicate that the measurement of cofactor precursor transport rates can inform the internal state of the cell when the external demands are altered or a causal genetic variation is introduced. Finally, genetic mutations that alter the clinical phenotype may also disrupt the 'natural' time scale hierarchy of interactions in the network

Systems analysis of host-parasite interactions.

Parasitic diseases caused by protozoan pathogens lead to hundreds of thousands of deaths per year in addition to substantial suffering and socioeconomic decline for millions of people worldwide. The lack of effective vaccines coupled with the widespread emergence of drug-resistant parasites necessitates that the research community take an active role in understanding host-parasite infection biology in order to develop improved therapeutics. Recent advances in next-generation sequencing and the rapid development of publicly accessible genomic databases for many human pathogens have facilitated the application of systems biology to the study of host-parasite interactions. Over the past decade, these technologies have led to the discovery of many important biological processes governing parasitic disease. The integration and interpretation of high-throughput -omic data will undoubtedly generate extraordinary insight into host-parasite interaction networks essential to navigate the intricacies of these complex systems. As systems analysis continues to build the foundation for our understanding of host-parasite biology, this will provide the framework necessary to drive drug discovery research forward and accelerate the development of new antiparasitic therapies

Systems biology of SNPs

Genome-scale networks can now be reconstructed based on high-throughput data sets. Mathematical analyses of these networks are used to compute their candidate functional or phenotypic states. Analysis of functional states of networks shows that the activity of biochemical reactions can be highly correlated in physiological states, forming so-called co-sets representing functional modules of the network. Thus, detrimental sequence defects in any one of the genes encoding members of a co-set can result in similar phenotypic consequences. Here we show that causal single nucleotide polymorphisms in genes encoding mitochondrial components can be classified and correlated using co-sets

Transcriptome-based reconstructions from the murine knockout suggest involvement of the urate transporter, URAT1 (slc22a12), in novel metabolic pathways.

URAT1 (slc22a12) was identified as the transporter responsible for renal reabsorption of the medically important compound, uric acid. However, subsequent studies have indicated that other transporters make contributions to this process, and that URAT1 transports other organic anions besides urate (including several in common with the closely related multi-specific renal organic anion transporters, OAT1 (slc22a6) and OAT3 (slc22a8)). These findings raise the possibility that urate transport is not the sole physiological function of URAT1. We previously characterized mice null for the murine ortholog of URAT1 (mURAT1; previously cloned as RST), finding a relatively modest decrement in urate reabsorptive capacity. Nevertheless, there were shifts in the plasma and urinary concentrations of multiple small molecules, suggesting significant metabolic changes in the knockouts. Although these molecules remain unidentified, here we have computationally delineated the biochemical networks consistent with transcriptomic data from the null mice. These analyses suggest alterations in the handling of not only urate but also other putative URAT1 substrates comprising intermediates in nucleotide, carbohydrate, and steroid metabolism. Moreover, the analyses indicate changes in multiple other pathways, including those relating to the metabolism of glycosaminoglycans, methionine, and coenzyme A, possibly reflecting downstream effects of URAT1 loss. Taken together with the available substrate and metabolomic data for the other OATs, our findings suggest that the transport and biochemical functions of URAT1 overlap those of OAT1 and OAT3, and could contribute to our understanding of the relationship between uric acid and the various metabolic disorders to which it has been linked

Network reconstruction of platelet metabolism identifies metabolic signature for aspirin resistance.

Recently there has not been a systematic, objective assessment of the metabolic capabilities of the human platelet. A manually curated, functionally tested, and validated biochemical reaction network of platelet metabolism, iAT-PLT-636, was reconstructed using 33 proteomic datasets and 354 literature references. The network contains enzymes mapping to 403 diseases and 231 FDA approved drugs, alluding to an expansive scope of biochemical transformations that may affect or be affected by disease processes in multiple organ systems. The effect of aspirin (ASA) resistance on platelet metabolism was evaluated using constraint-based modeling, which revealed a redirection of glycolytic, fatty acid, and nucleotide metabolism reaction fluxes in order to accommodate eicosanoid synthesis and reactive oxygen species stress. These results were confirmed with independent proteomic data. The construction and availability of iAT-PLT-636 should stimulate further data-driven, systems analysis of platelet metabolism towards the understanding of pathophysiological conditions including, but not strictly limited to, coagulopathies

Quality of Radiomic Features in Glioblastoma Multiforme: Impact of Semi-Automated Tumor Segmentation Software.

ObjectiveThe purpose of this study was to evaluate the reliability and quality of radiomic features in glioblastoma multiforme (GBM) derived from tumor volumes obtained with semi-automated tumor segmentation software.Materials and methodsMR images of 45 GBM patients (29 males, 16 females) were downloaded from The Cancer Imaging Archive, in which post-contrast T1-weighted imaging and fluid-attenuated inversion recovery MR sequences were used. Two raters independently segmented the tumors using two semi-automated segmentation tools (TumorPrism3D and 3D Slicer). Regions of interest corresponding to contrast-enhancing lesion, necrotic portions, and non-enhancing T2 high signal intensity component were segmented for each tumor. A total of 180 imaging features were extracted, and their quality was evaluated in terms of stability, normalized dynamic range (NDR), and redundancy, using intra-class correlation coefficients, cluster consensus, and Rand Statistic.ResultsOur study results showed that most of the radiomic features in GBM were highly stable. Over 90% of 180 features showed good stability (intra-class correlation coefficient [ICC] ≥ 0.8), whereas only 7 features were of poor stability (ICC < 0.5). Most first order statistics and morphometric features showed moderate-to-high NDR (4 > NDR ≥1), while above 35% of the texture features showed poor NDR (< 1). Features were shown to cluster into only 5 groups, indicating that they were highly redundant.ConclusionThe use of semi-automated software tools provided sufficiently reliable tumor segmentation and feature stability; thus helping to overcome the inherent inter-rater and intra-rater variability of user intervention. However, certain aspects of feature quality, including NDR and redundancy, need to be assessed for determination of representative signature features before further development of radiomics

iAB-RBC-283: A proteomically derived knowledge-base of erythrocyte metabolism that can be used to simulate its physiological and patho-physiological states

<p>Abstract</p> <p>Background</p> <p>The development of high-throughput technologies capable of whole cell measurements of genes, proteins, and metabolites has led to the emergence of systems biology. Integrated analysis of the resulting omic data sets has proved to be hard to achieve. Metabolic network reconstructions enable complex relationships amongst molecular components to be represented formally in a biologically relevant manner while respecting physical constraints. <it>In silico </it>models derived from such reconstructions can then be queried or interrogated through mathematical simulations. Proteomic profiling studies of the mature human erythrocyte have shown more proteins present related to metabolic function than previously thought; however the significance and the causal consequences of these findings have not been explored.</p> <p>Results</p> <p>Erythrocyte proteomic data was used to reconstruct the most expansive description of erythrocyte metabolism to date, following extensive manual curation, assessment of the literature, and functional testing. The reconstruction contains 281 enzymes representing functions from glycolysis to cofactor and amino acid metabolism. Such a comprehensive view of erythrocyte metabolism implicates the erythrocyte as a potential biomarker for different diseases as well as a 'cell-based' drug-screening tool. The analysis shows that 94 erythrocyte enzymes are implicated in morbid single nucleotide polymorphisms, representing 142 pathologies. In addition, over 230 FDA-approved and experimental pharmaceuticals have enzymatic targets in the erythrocyte.</p> <p>Conclusion</p> <p>The advancement of proteomic technologies and increased generation of high-throughput proteomic data have created the need for a means to analyze these data in a coherent manner. Network reconstructions provide a systematic means to integrate and analyze proteomic data in a biologically meaning manner. Analysis of the red cell proteome has revealed an unexpected level of complexity in the functional capabilities of human erythrocyte metabolism.</p

Systems biology analysis of drivers underlying hallmarks of cancer cell metabolism.

Malignant transformation is often accompanied by significant metabolic changes. To identify drivers underlying these changes, we calculated metabolic flux states for the NCI60 cell line collection and correlated the variance between metabolic states of these lines with their other properties. The analysis revealed a remarkably consistent structure underlying high flux metabolism. The three primary uptake pathways, glucose, glutamine and serine, are each characterized by three features: (1) metabolite uptake sufficient for the stoichiometric requirement to sustain observed growth, (2) overflow metabolism, which scales with excess nutrient uptake over the basal growth requirement, and (3) redox production, which also scales with nutrient uptake but greatly exceeds the requirement for growth. We discovered that resistance to chemotherapeutic drugs in these lines broadly correlates with the amount of glucose uptake. These results support an interpretation of the Warburg effect and glutamine addiction as features of a growth state that provides resistance to metabolic stress through excess redox and energy production. Furthermore, overflow metabolism observed may indicate that mitochondrial catabolic capacity is a key constraint setting an upper limit on the rate of cofactor production possible. These results provide a greater context within which the metabolic alterations in cancer can be understood

A detailed genome-wide reconstruction of mouse metabolism based on human Recon 1

<p>Abstract</p> <p>Background</p> <p>Well-curated and validated network reconstructions are extremely valuable tools in systems biology. Detailed metabolic reconstructions of mammals have recently emerged, including human reconstructions. They raise the question if the various successful applications of microbial reconstructions can be replicated in complex organisms.</p> <p>Results</p> <p>We mapped the published, detailed reconstruction of human metabolism (Recon 1) to other mammals. By searching for genes homologous to Recon 1 genes within mammalian genomes, we were able to create draft metabolic reconstructions of five mammals, including the mouse. Each draft reconstruction was created in compartmentalized and non-compartmentalized version via two different approaches. Using gap-filling algorithms, we were able to produce all cellular components with three out of four versions of the mouse metabolic reconstruction. We finalized a functional model by iterative testing until it passed a predefined set of 260 validation tests. The reconstruction is the largest, most comprehensive mouse reconstruction to-date, accounting for 1,415 genes coding for 2,212 gene-associated reactions and 1,514 non-gene-associated reactions.</p> <p>We tested the mouse model for phenotype prediction capabilities. The majority of predicted essential genes were also essential in vivo. However, our non-tissue specific model was unable to predict gene essentiality for many of the metabolic genes shown to be essential in vivo. Our knockout simulation of the lipoprotein lipase gene correlated well with experimental results, suggesting that softer phenotypes can also be simulated.</p> <p>Conclusions</p> <p>We have created a high-quality mouse genome-scale metabolic reconstruction, iMM1415 (<it>Mus Musculus</it>, 1415 genes). We demonstrate that the mouse model can be used to perform phenotype simulations, similar to models of microbe metabolism. Since the mouse is an important experimental organism, this model should become an essential tool for studying metabolic phenotypes in mice, including outcomes from drug screening.</p